Crossover from Kramers to phase-diffusion switching in hysteretic DC-SQUIDs 
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We have measured and propose a model for switching rates in hysteretic DC-SQUID in the regime 
where phase diffusion processes start to occur. We show that the switching rates in this regime are 
smaller than the rates given by Kramers' formula due to retrapping of Josephson phase. The 
retrapping process, which is affected by the frequency dependent impedance of the environment of 
the DC-SQUID, leads to a peaked second moment of the switching distribution as a function of 
temperature. The temperature where the peaks occur are proportional to the critical current of the 
DC- SQUID. 

PACS numbers: 74.40.+k, 85.25.Dq, 74.50 +r 



The phenomenon of phase-diffusion in Josephson junc- 
tions has been extensively studied over the past 20 years 
0, 0,111 0]- These studies, in general, focused on deep 
sub-micron junctions with very low critical currents, I c , 
so that the Josephson coupling energy Ej was of the or- 
der fcflT . Extensive diffusion of the phase typically oc- 
curred before the junction switched to the so-called run- 
ning state with voltage near the superconducting gap. 
With the advent of research on flux qubits for quan- 
tum computing using SQUID magnetometers with much 
larger, unshunted junctions for readouts (J c > 1/iA), 
this phase-diffusion process has reappeared, rather dra- 
matically, in a different regime. In these larger junc- 
tions, which are the focus of this work, the crossover 
from simple Kramers activation to phase-diffusion be- 
fore switching tends to occur at higher temperatures and 
involves only a small number of phase slip events be- 
fore the switching process. Nevertheless the effects of 
this crossover can be quite pronounced, appearing, per- 
haps most dramatically, as a significant narrowing of the 
switching distribution and thus, in a sense, a greatly in- 
creased sensitivity of the magnetometer, albeit, at the 
expense of greater back action on the qubit. 

The standard picture for the dynamics of a Josephson 
junction, shown in Fig. la, is that of a particle with a 
mass proportional to the junction capacitance moving in 
a potential U(tp) = —Ejcos(tp) — It,ip, where ip is the 
phase and lb is the bias current of the junction. For suf- 
ficiently low damping around the plasma frequency (oj p ) 
of the Josephson junction, the energy gain as the parti- 
cle moves from one barrier to the next will exceed the 
loss due to damping. So, switching from the supercur- 
rent state to running state is triggered by a single event, 
the phase particle escaping from a potential well, and no 
phase diffusion occurs. This situation is usually realized 
in large area Josephson junctions. On the other hand, 
if the damping around the junction plasma frequency is 
sufficiently high then at low bias the phase particle will 
always retrap in a local minima after escape and a finite 
resistance phase-diffusion branch exists on the I-V curve 



of the junction. This case is usually realized in ultra- 
small Josephson junctions and can be observed down to 
lowest temperatures attainable with a dilution refrigera- 
tor ■ For low damping, the switching rate equals to the 
thermally activated (TA) escape rate from the potential 
minima which is given by well known Kramers' formula 
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where fi is the effective attempt frequency (f2 « u p ), AU 
is the height of potential barrier from a local minimum to 
the next maximum (see also Fig. la) and a t is a damping 
dependent factor (0 < a* < 1) 0. For the case of high 
damping, where phase diffusion is always present at low 
bias, a similar expression can be given for the switch- 
ing rate although the meaning of the potential barrier 
and attempt frequency are very different 0. For inter- 
mediate damping a third regime occurs that bridges the 
previous two. At lower temperatures a single escape can 
trigger the switching of the junction from the metastable 
minima but at higher temperatures, where the mean es- 
cape current is reduced, the energy gain is less and 
phase diffusion occurs. In this case the switching rate is 
not only determined by the escape rate but also by the 
probability of retrapping of the phase particle after it has 
escaped. 

The total transition rate from metastable minima 
to running state can be written as a sum of rates 
through different escape and retrapping sequences as 
r = T,n=oP {n)T{n) - Here r(n) stands for the switching 
rate of a process where the particle retraps n times before 
the runaway starts and p^ for the probability of such a 
process, p( n > can be expressed in terms of the retrapping 
probability (P r t), which we define as the probability that 
after the ascent of the phase particle to a local maximum 
it will be retrapped in one of the subsequent minima 
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For the case relevant to our 
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FIG. 1: a) Schematics of washboard potential of DC-SQUID. 
Angular frequency of small oscillations of the fictitious phase 
particle is ui p . To escape phase particle has to overcome po- 
tential barrier AU. b) Schematics of circuit used in calcula- 
tion of retrapping probability. R s , Cb model high frequency 
damping seen by a Josephson junction with critical current 
I c . I„ represents noise and h bias currents, respectively. C 
marks capacitance of the junction, c) In experiments single 
junction is replaced by a DC-SQUID which effective critical 
current I c (&b) can be modulated by external magnetic field. 

data, where Tt <C oj p , this gives 

r = r T (i - p rt ) E ^ = r T (i - p rt ) ln(1 "/ rt) " 1 

71= 

(2) 

In Eq. 2, the switching rate factorizes into two parts, one 
of which is the Kramers rate (1) and the other a function 
of the retrapping probability. Since the latter is always 
less than or equal to one, the actual switching rate is 
smaller than or equal to the rate predicted by Kramers 
equation (1). At zero temperature, consideration of the 
balance between energy gain from the bias source and 
energy losses due to the damping leads to the conclusion 
that P r t(h) will abruptly switch from 1 to at some bias 
current I r — the retrapping current. For ohmic damping 
and Q > 1 (Q = uo p CR), I rQ = 4I c /irQ 00. 

To find P r t for finite temperatures and frequency de- 
pendent damping, we rely on Monte Carlo simulations. 
We calculate P rt for the RCSJ model with a series RC 
shunt added to account for the stronger high frequency 
damping (Fig. l.b), e.g. due to the leads. Extensive dis- 
cussion of this model along with equations of motion for 
the phase can be found in ref. 0- The damping at to p is 
assumed to be much stronger than that at low frequency 
as required for the coexistence of phase-diffusion and hys- 
teresis. However, the exact value of the low frequency 
damping is not important here, so we omit the quasipar- 
ticle resistance and bias resistors which are much larger 
than the typical transmission line impedance (~ 100 f2) 
seen by the junction near lo p . To obtain P rt , the parti- 
cle is initiated repeatedly at the top of a barrier and its 
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FIG. 2: Retrapping probability as a function of normalized 
bias current for circuit shown in Fig. 1. Parameters used 
in calculations are I c = 2.90 /LtA, C = 260 fF, R s = 75 fl, 
Ct = 15 pF corresponding to sample B. Curves from left to 
right correspond to temperatures (dashed line), 0.1, 0.5, 1.0, 
2.0, 5.0 K. Inset shows retrapping probability vs. normalized 
bias current for R s = 125, 100, 75, 50 il (from left to right) 
at T = 2.0 K. I c , C, Cb have been chosen the same as in the 
main Figure. 



subsequent motion in the presence of Nyquist noise is 
monitored until it is retrapped in one of the metastable 
minima or switches to the running state. 

Figure 2 shows this calculated P r t as a function of bias 
current for a range of temperatures that is important for 
our data. At zero temperature, P rt switches abruptly, as 
expected from energy balance. At higher temperatures, 
however, the retrapping probability curve broadens and 
shifts to higher bias currents with i|0% x yi/2 wnere 
7 b 50% is defined as P rt (/ fc 50% ) = 0.5. Even though the 
overall shift is toward higher bias as T increases, P r t(Ib) 
also decreases slightly from unity for If, < I r where I r is 
the actual retrapping current for the circuit model in Fig. 
la, i.e. I r = I r (R s ,Cb) calculated at T — 0. For I\, < I r , 
the simulations show 1 — P rt decreases essentially expo- 
nentially as lb decreases. For this to occur, the frequency 
dependence of the damping is crucial. If thermal fluctu- 
ations give enough initial velocity to the particle, the 
linear component of -J£ will increase, lowering the veloc- 
ity dependent damping and thereby permitting switching 
even though lb < I r - Due to high probability of retrap- 
ping the switching rate T(Ib) deviates substantially from 
Kramers' rate Tt(/&) in this range of bias currents. 

The effects of the choices of the circuit parameters Cb 
and R s have been investigated using the simulations. The 
inset in Fig. 2 shows P r t(Ib) for a range of R s . As 
can be seen, P r t{Ib) shifts to lower bias currents with 
j50% x \Jji s as predicted by I r w AI c /ttQ p where for our 
non-Ohmic circuit model Q p = lo p CR s . We have chosen 
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TABLE I: Total critical current (Ico), sum of capacitances of 
two junctions (C), inductance (L) and assymmetry of critical 
currents of the two junctions of the DC-SQUID (a). 



DC-SQUID 


ho (a*A) 


C (fF) 


L (pH) 


a 


A 


4.25 


90 


37 


0.025 


B 


2.90 


260 


70 


0.05 



values of Cb to be sufficiently large (15-20 pF) that the 
damping is nearly frequency independent for frequencies 
near ui p . For these and larger values of Cb, P r t(h) is 
essentially independent of Cb for Jfe > I r . Only in the re- 
gion lb < I r , does its exact value become important. For 
determining the onset of phase-diffusion as T increases, 
the exact value of Cb is not critical. However, the details 
of the regime where more extensive phase diffusion takes 
place do depend on the exact frequency dependence of 
the damping. 

The samples we used in this study were DC-SQUIDs 
fabricated using a self-aligned Nb trilayer process [Icj . 
The inductances of the SQUIDs were chosen small 
enough that (3l = 2ttLI c q/^q << 1. So, for most flux 
biases, the potential of the SQUIDs could be approx- 
imated well by that of a single junction with an effec- 
tive critical current I c ($b) (Fig. lc). Two samples with 
different parameters were measured in our two laborato- 
ries. The results of these two measurements are essen- 
tially identical. The parameters of the samples are listed 
in Table 1. We determined the total critical current of 
the SQUIDs by fitting the switching rates in the tem- 
perature interval 0.4-1.0 K, where, for <f>b = 0, the rates 
followed the Kramers rate (Eq. 1) well. We estimate the 
total capacitance (C) of two junctions of the SQUID from 
their area and the specific capacitance of 45 fF /Vm 2 and 
calculate the inductances (L) by using a 3D inductance 
calculation program (llj . 

The sample is placed in a copper can filled with He ex- 
change gas and located on a temperature regulated stage 
of a dilution refrigerator. Extensive electrical filtering 
is used to prevent interference and higher temperature 
thermal noise from reaching the sample. The bias cur- 
rent of the SQUID is ramped at constant rate. When a 
voltage appears across the device, the value of bias cur- 
rent at which switching occurred is recorded. This is 
typically repeated over 10 4 times at a given temperature 
and flux bias and a histogram of the switching currents 
is compiled. In Fig. 3, typical switching rates versus the 
ratio of barrier height to temperature from sample A are 
plotted. Here the barrier height is calculated from the 
full two dimensional (2D) potential of the DC-SQUID 
for each value of bias current. According to Eq. (1) 
the rates should fall onto one straight line in this plot 
if the ln(wp) dependence of rate on the bias current is 
neglected. In the temperature range 0.4. .1.8 K this is in- 
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FIG. 3: Switching rate of DC-SQUID A as a function of bar- 
rier height to temperature ratio (solid lines) for temperatures 
0.4, 1.8, 2,2, 2.6, 3.0 K and flux bias of DC-SQUID $& = 0.167 
$o. Dash-dot line shows prediction of Eq. (1). Dashed lines 
represent prediction of Eq. (2) for temperatures 1.8, 2.2, 2.6, 
3.0 K. Note that predictions by eq. (1) and (2) coincide for 
temperature 1.8 K and lower. Retrapping probability is calcu- 
lated for the model with parameters R s — 70 SI and Cb ~ 20 
pF. Inset shows measured switching rate at temperature 2.2 
K (solid line). Dashed line show calculated rates for R s = 60, 
70 and 80 SI. 

deed the case. At higher temperatures and larger value of 
AU/ksT, however, the switching rates increasingly devi- 
ate from the simple exponential dependence described by 
Eq. (1). Following our earlier argument we expect this to 
happen when retrapping takes place. To verify this, we 
also plot the modified rate calculated from simulations 
of Eq. (2) in the same figure. To achieve a reasonable 
speed, a ID approximation to the full 2D potential with 
an effective I c ($b) is used to obtain P rt . The agreement 
between measurement and the model is very good despite 
the simplicity of the circuit model. 

Next we study the crossover to phase diffusion as a 
function of temperature and / c ($b). For clarity, in this 
large data set, only the first two moments of switching 
histograms, lb and a], rather than the complete T(Ib) 
curves, are presented. Figure 4. shows lb and <jj for 
different values of I c ($b) f° r sample B. One can notice 
three distinct regions in the temperature dependence of 
the width 07. For flux bias $b « and T < 80 mK the 
escape is by quantum tunneling (QT) giving a width that 
is nearly independent of T. For higher T , in the range of 
0.1 — 2.0 K, and for larger values of I c (&b) 07 agrees well 
with the predictions of TA (Eq. l).The cross-over tem- 
perature to QT decreases with J c ($f,), as expected, with 
the very lowest I c data never completely reaching the QT 
regime. At still higher temperatures, we observe a third 
regime where the width of the distribution decreases as 
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FIG. 4: Mean (top panel) and width (bottom panel) of 
switching distribution as a function of temperature for sample 
B. Solid diamonds correspond to flux bias <!>(, = $o, solid 
rectangles to $ 6 = 0.239 <E>o, solid triangles to $ b = 0.324 $o, 
empty diamonds to $6 = 0.409 $o and empty rectangles to 
$(, = 0.474 "I>o- Dashed lines show mean and width calculated 
using rate given by Eq. (1). Solid lines represent results of 
calculation based on Eq. (2). Inset shows temperature where 
the width of switching distribution peaks as a function of ef- 
fective critical current of the DC-SQUID (circles) and results 
of calculation based on Eq. (2) (triangles). Ramp rate of bias 
current is 88 fiA/s. 

temperature is increased . Based on our earlier discus- 
sion, this results from retrapping of the system after its 
initial escape 0. Similarly to the second moment, lb 
agrees well with theory in the QT and TA regimes. It is 
greater than predicted by TA after the onset of retrap- 
ping since a greater tilt to the potential is required to 
overcome the retrapping and initiate switching. 

Figure 4 also shows calculations based on Eq. (2) (solid 
lines). As one can see Eq. (2) predicts well the anomalous 
temperature dependence of the width for <£>/, = 0, 0.239 
and 0.324 <& but gives only a fair fit for $ b = 0.409 
and 0.474 <5> . For $ h « 0.5 $ , the ID potential of a 
single junction with effective critical current I c (<&&) is a 
rather poor approximation to the actual 2D potential of 



the DC-SQUID which can explain the discrepancy. Nev- 
ertheless, for all values of effective critical currents, cal- 
culations based on Eq. (2) predict reasonably accurately 
the temperature T*at which the width of the switching 
distribution is a maximum. The inset of Fig. 4. shows T* 
as a function of the effective critical current along with 
predictions of Eq. (2). As this figure shows, T* (J c ) spans 
more than one order of magnitude for a given SQUID and 
scales linearly with the effective critical current. Our cal- 
culations also show that the observed linear T* vs. I c ($b) 
dependence follows only for sufficiently small values of 
high frequncy damping, i.e. large R s . For larger damp- 
ing, this dependence deviates from linear and approaches 
quadratic. The simulations of dynamics of the phase par- 
ticle (near <!>(, = 0) show also that at the temperature 
T* the retrapping probability around the mean of the 
switching distribution is relatively low (P r t(h) ~ 0.1), 
so most escapes of the phase particle from local min- 
ima lead to switching. In those cases when retrapping 
occurs, the number of wells the phase particle moves be- 
fore retrapping is, on average 2.. 3. As the temperature 
rises further and the average switching current decreases 
advances of the phase particle before retrapping become 
shorter due to smaller energy input from bias source but 
more frequent due to the increased level of current fluc- 
tuations. This leads eventually to the appearance of a 
stable phase diffusion state where the phase particle dif- 
fuses down the tilted Josephson potential with uniform 
average speed giving rise to a measurable phase diffusion 
voltage. We observe such a phase diffusion voltage (~1 
fiY) for $ b = 0.409 and 0.474 <J> (sample B) at T « 4.2 
K. 

In summary, we observe the crossover from Kramers 
escape to escape affected by phase diffusion in relatively 
large Josephson junctions having I c of the order of 1 fxA. 
This is manifested by a peak in the width of the switching 
distribution 07 vs T at a temperature T*(J C ) which scales 
linearly with I c . For T > T*(I C ), ai can decrease signif- 
icantly below that expected for Kramers escape. Our 
data are in good agreement with the results of Monte 
Carlo simulations based on a simplified circuit model for 
the junction in which the higher damping at the plasma 
frequency is represented by a series RC shunt. 
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